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ABSTRACT 

Using the IRAM 30m telescope, we have detected the '^CO 7 = 2 - 1, 4 - 3, 5 - 4, and 6-5 emission lines in the millimeter- 
bright, blank-field selected AGN COSMOS J100038+020822 at redshift z = 1.8275. The sub-local thermodynamic equilibrium 
(UTE) excitation of the J = 4 level implies that the gas is less excited than that in typical nearby starburst galaxies such as NGC253, 
and in the high-redshift quasars studied to date, such as Jl 148-1-5251 or BR1202-0725. Large velocity gradient (LVG) modeling of the 
CO line spectral energy distribution (CO SED; flux density vs. rotational quantum number) yields Ht densities in the range 10''^ - lO*" 
cm"^, and kinetic temperatures between 50 K and 200 K. The H2 mass of (3.6-5.4)x 10'" Mq implied by the line intensities compares 
well with our estimate of the dynamical mass within the inner 1.5 kpc of the object. Fitting a two-component gray body spectrum, we 
find a dust mass of 1.2 x 10' Mq, and cold and hot dust temperatures of 42+5 K and 160+25 K, respectively. The broad Mgll line 
allows us to estimate the mass of the central black hole as 1.7 x 10' Mq. Although the optical spectrum and multi- wavelength SED 
matches those of an average QSO, the molecular gas content and dust properties resemble those of known submillimeter galaxies 
(SMGs). The optical morphology of this source shows tidal tails that suggest a recent interaction or merger. Since it shares properties 
of both starburst and AGN, this object appears to be in a transition from a strongly starforming submillimeter galaxy to a QSO. 

Key words. Galaxies: evolution - Galaxies: starburst - Galaxies individual (J100038+020822) - Galaxies: ISM - Galaxies: quasars: 
emission lines - Galaxies: high-redshift 
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1 . Introduction serva tions dSanders et al.lll988at iPage et alj|2004l: IStevens et all 

> • I2005h . 

Submilhmeter blank field surveys have discovered a popu- The u biquitous presence of AGN activity in most SMGs 

latign of dust enshrouded high-redshiftg^axies d SmailetalJ ^Alexand er et al. 2005) suggests a close link between the AGN 

[19971 [ Hughes etal.| [1998|; [Bales etal.| [1999D, which are mas- starburst activity at redshifts z ~ 1 - 3. If a QSO is found to 

sive systems with huge molecular gas reservoirs andstar f^^.j^ i^^i^inous, it means that large amounts of gas and dust 

formation at high rates Neri etal.[ [20031 i GreveetalJ [20051 should still be present. Such a source may, in fact, be a good can- 

[Solomon & Vanden BoutI l2005[). Their faint X-ray emission ^^^^^^ ^„ ^^-^^^ i„ transition from a starburst to a QSO, 

([Alexander et al. 2005) suggests that most of the submillimeter jn particular when it shows absorbed X-ray emission ( PagilFd] 

output is not powered by active galactic nuclei (AGN), but by ^ IStevens et al.l 12001 . Yet only a few high-redshift com- 

massive star formation. These starbursting submillimeter galax- ^^^^^^ starburst/AGN have been studi ed in its molecula r and 

les (SMGs) account for a substancial fraction of the far-infrared muitj-waveleng th con tinuum emission (lRowan-Robinsonl [2000t 

(IR) background, and cuiTent models suggest this population |Le Floc'h et alJ[205% ICoppin et al.ll2008l) . 

may represent the formation of massive spheroidals at high- ^, * u /r^m^ ■ 1 ■ 

redshift do nl Tj [ 200lh Observations of carbon monoxide (CO) in galaxies are im- 

" ^ portant probes of the physical conditions of the cold and warm 

In the local Universe, an evolutionary connection be- molecular gas in the galactic nuclei and disks. They provide esti- 

tween starbursting ultralummous infrared galaxies (ULIRGs) ^ates of the total amount of gas available to fuel starburst and/or 

and QSOs has been suggested ( Sanders et al. 19883) and dis- ^GN activity, and the CO line profile and intensity can be used to 

cussed controversially (e.g. Sanders et al. 1988a b; GenzeleLaU obtain important information about the galaxy ki nematics, such 

[1998t [ Tacconi etalj |2002). This evolutionary cycle has found as dynamical mass or size of the emitting region dSolomon et al.1 

support from hydrodvnamical simulations of galaxy formation |T997[ [Solomon & Vanden Bo^[2005[) . 

(e.g. [ DiMatteoetalJ [ 2005 fc [ Hopkins et al. [[ 2005i [2006|) and ob- g^^^^^^ jl^^i^ diagnostic value, great efforts have been 

made to observe CO emission lines in SMGs. These stud- 

Send offprint requests to: M. Aravena ies have benefited from deep radio continuum imaging (e.g. 
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VLA 1.4 GHz) to locate the SMG accurately dlvison et alj 
I2OO2I I2005L I2OO7 I), and from the determination of optical 
spectroscopic redshifts (e.g. Chapman etal. 2005). Due the 
small spectroscopic bandwidths of current (sub)millimeter tele- 
scopes and interferometers, this was necessary to permit the 
proper frequency tuning for line observations (e.g. of CO, 
rCIl rem). Y et only 19 SMGs have been reported in CO so 
far (iFraver e t al. 1998, 1999; Andreani et al. 2000; Sheth et al) 
20041: iHa inline et al. 2004; Neri et all 120031; [C reve et al. 2005; 
Frayer et al.,.2008. : .Coppin et al...2008t iTacconi et aL,200 8) and 



only a few have been observed in multiple molecular tran- 
sitions to stud y the excitation conditions of their molecular 
gas reservoir (I Solomon & Vanden Boud l2005h . Spatial struc- 
ture and dynamics were studied for some SMGs th rough high 
resolution CO imaging (iTacconi et all l2006l l2008i) . Overall, 
only ~ 50 high-redshift (z > 1) objects have been detected 
in CO, most of wh i ch are luminous optically se l ected QSOs 
dOrnont et all 119961; iGuilloteau et al.l Il997lll999l; ICariUi et al.1 



2002HWalter et al.l2003HBertoldi et al.l200 3t lBeelen et al.l2004l: 



Riechers et all 120061: JCarilli et al.' '2007') and high-reds hift ra- 
dio galaxies (HzRG; |^ Breuck et al. 2003, 2005; Klam er et all 
I2005t iPapadopoulos eral.li2005h . 

Here we report the detection of CO 2 - 1,4-3, 5 - 4, and 
6-5 line emission from a millimeter selected QSO, J100038.01 + 
020822.4. Its CO line intensities and flux density ratios allow 
us to estimate its molecular gas content and its excitation con- 
ditions. Its spectral properties suggest that this object may be 
evolving from a starburst to a QSO. 

After a description of the observations in Section 2, we 
present the molecular gas and dust properties of our source 
in Sections 13.11 to 13.41 We study its morphology and multi- 
wavelength properties in Sections 13.51 13.61 and 13.71 We discuss 
the results in Section 4 and give a brief summary in Section 
5. Throughout, we use a ACDM cosmology, Hq - 70 km s"' 
Mpc ', Qa = 0.7 and Qm = 0.3. 



2. Observations 

2.1. Source selection 

Pan-chromatic surveys are essential for an understanding of 
galaxy properties and their evolution. The cosmic evolution sur- 
vey (COSMOS) is the first multi-wavelength survey to cover 
a sufficiently large area (1.4° x 1.4°) at appropriate depth over 
nearly the entire electromagnetic spectrum to provide a compre- 
hensive view of gal axy formation and large scale structure (see 
IScoviUe et al.ll2007h . 

As part of the COSMOS project, Bert oldi et all (|2007|) 
mapped the central ~ 20' x 20' of the COSMOS field at 1.2 mm 
(250 GHz) using the Max-Planck millimeter bolometer cam- 
era (MAMBO) at the IRAM 30m telescope. An accompany- 
ing VLA radio imagi ng project (1.4 GHz to ~ 10 juJy rms; 
ISchinnerer et alj 12007) allowed for the identification of radio 
counterparts for 24 millimeter sources. 

The MAMBO source J100038.01-t-020822.4 (hereafter: 
J100038) is one of the strongest millimeter selected sources in 
COSMOS (boosting corrected S 1.2 mm - 4.6 ± 0.9 mJy). 

It hosts an X-ray luminous (Lx ~ lO"*"* erg s"'), absorbed 
(logA^(H) ~ 22 - 23 cm'^) AGN that classifies it as an ob 



scured QSO (iBrusa et alJl2007l;lMainieri et alj|2007[). and its op- 
tical s pectrum is typical of a broad line (BL) AGN (^ Trump et alJ 
I2007h . The relatively faint 1.4 GHz radio emission (5i.4ghz — 
237 + 27 fj.}y) suggests that the millimeter emission arises from 
a starburst. In fact, the radio-to-millimetre flux ratio can be used 



as a redshift indicator (ICarilli&Yunll 19991) that implies z = 1.9, 
which is consistent with the spectroscopic redshift. Since this 
i s a millimeter sele cted QSO that has been detected in X-rays 
dStevens et alJ2005l) . it likely constitutes a transitional case from 
SMG to QSO. 

Three independent spectroscopic re dshift measu r ement s 
exist for J10003 8: 1.8325 + 0.0023 (Trump et al.' '2007|), 
1.825 + 0.002 (P rescott et al. 2006) and 1.8289 + 0.002 (Marco 
Scodeggio, priv. comm.), which we have averaged to z = 
1.8288 + 0.0037. 

2.2. CO observations 

The CO observations were performed with the IRAM 30m tele- 
scope on Pico Veleta, Spain, during the winter 2006/2007 in 
good millimeter weather conditions (precipitable water vapor 
< 5 mm). We observed the CO 2 - 1 (redshifted to 81.551 GHz) 
and CO 5 - 4 (203.850 GHz) lines simultaneously using the A/B 
receiver configuration, and the CO 4 - 3 (163.088 GHz) and CO 
6-5 (244.603 GHz) lines using the C/D receiver configuration. 

The observations were centered on the optical Subaru / band 
position of J100038 at lO'' 00"' 38.01' -h02° 08' 22.6" (J2000). 
The optical position lies 0.04" away from the VLA 1.4 GHz 
position. The beam size of the IRAM 30m telescope at 200 GHz 
and 81 GHz is 12" and 30", respectively. 

We observed in wobbler switching mode with a wobbler rate 
of 0.5 Hz and a wobbler throw of 60" in azimuth. Due to its 
proximity to the source, we used Saturn as main pointing and 
focus calibrator The pointing was checked every hour and was 
stable within 3" during all runs. 

Typical values for the system temperatures were ~ 140 K, 
400 K, 460 K and 600 K for the 3 mm, 2 mm, and lower and 
higher 1 mm bands, respectively. We calibrated every 12 minutes 
with hot/cold absorbers and estimate the fluxes to be accurate to 
+ 10% at 3 and 2 mm and 20% at 1 mm. As spectrometers we 
used the 512 x 1 MHz filterbanks for the 3 mm receivers, and 
the 256 X 4 MHz filterbanks for the 2 and 1 mm receivers (1 
GHz bandwidth). We reduced the data using CLASS, removing 
scans with strongly distorted baselines, substracted linear base- 
lines in the remaining scans and rebinned the averaged spectra 
to velocity resolutions of 100, 120, 90 and 130 km s"' for the 
CO 2 - 1,4-3, 5-4 and 6-5 lines, respectively. This led 
to baseline antenna temperature (T^) rms noise levels of 0.04, 
0.17, 0.44 and 0.4 mK. To convert to main beam temperatures, 
Tmb, we multiply by the ratio between the forward and beam 
efficiencies, Feif/Seff, at the observed frequencies. Flux densities 
were obtained using the conversion factor Sv/T,^\, = 4.95 Jy/K 
for the IRAM 30m telescopefl 

3. Results and Analysis 

Figure [T] shows the resulting spectra and Table [1] summarizes 
the measured line parameters. The line profiles for the different 
transitions appear Gaussian and similar to each other, with no 
signs of velocity structure, as might be expected for a merger an d 
is often found for SMGs dGreve et al.ll200l IWeiB et al.llBo5al) . 
However, the signal to noise ratio is too low to examine the in- 
dividual line profiles. Although the CO 4-3 line appears to 
be slightly shifted in velocity relative to the other lines, the line 
widths are very similar, with an average of Avfwhm = 417 + 48 
km s The spectral energy distribution of the CO rotational 
emission (CO SED) peaks at the 7 = 5-4 transition (Fig. 2 and 



http : //www . iram . es/IRAMES/telescope . html 
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Fig. 1. Observed spectra of the CO 7 = 2 - 1, 4 - 3, 5 - 4 and 6 - 5 emission lines in J100038+020822, obtained with the IRAM 
30m telescope. Gaussian fits to the spectra are shown as dotted lines. The velocity reference is at z = 1.8275. 



Table 1. Observed line parameters for J100038. 



Transition 


Vobs " 


HPBW " 




'/ " cl 
J nib 




vfwhm 


Ico' 


^co 






[GHz] 


["] 


[mK] 


[mK] 


[mJy] 


[km s-'] 


[Jy km s"'] 


10'" [Kkms-' pc^] 


lO'* [Lq] 


C02- 1 


81.551 


30.2 


0.51 


0.62 


3.08 


406 ±31 


1.33 


5.6 ±0.35 


0.22 ± 0.01 


C04-3 


163.088 


15.1 


1.18 


1.69 


8.39 


427 ± 73 


3.81 


40 ±0.51 


1.25 ±0.16 


CO 5 -4 


203.850 


12.1 


1.46 


2.33 


11.56 


397 ± 161 


487 


3.3 ±0.8 


2.00 ± 0.49 


C06-5 


244.603 


10 


1.19 


2.09 


10.35 


443 ± 75 


486 


2.3 ± 0.4 


2.40 ± 0.47 



Notes: " Observed frequency; * Half power beam width of the IRAM 30m telescope at the observed frequency; Antenna temperature; Main 
beam temperature; Flux density; CO line full width half maximum; * Integrated CO intensity ; '' CO luminosity; ' Integrated CO luminosity. 



Table [U. Using the detected CO lines we estimate a systemic 
CO redshift z - 1.8275 + 0.0013, which we use throughout this 
paper. 

3.1. Large Velocity Gradient model 

To study the molecular gas excitation in J100038 we compared 
the relative CO line intensities with those predicted by a sin- 
gle component large velocity gradient (LVG) model, as suming 
sphe rical geometry. We used the collision rates from iFlowed 
(1200 li) with an ortho-para H2 ratio of 3 and a fixed CO abun- 
dance per velocity gradient of [CO ]/(dv/dr) - 10"^ pc (km 
s-^r' re.g. lWeiB et a'l]l2005bll2()07h . 

The best fit to the data is provided by a model with a kinetic 
temperature, Tkin of 95 K and an H2 density, n(H2) of 10^^ cm"^ 
(model A, Fig.|2]i. Similarly good fits are achieved with Tidn = 60 
K and n(H2) = 10^^ cm"^ (model B), or Tidn = 200 K and «(H2) 
= 10^-^ cm-^ (model C). They all show a peak in the CO SED at 
the 7 = 5-4 line, but predict different higher-/ intensities. 



Figure [3] shows the logj^^^ distribution as a function of the 
kinetic temperature and H2 density. The temperature-density pa- 
rameter space is degenerate at densities between lO-' " cm"-' and 
lO'* " cm"^. The most likely ranges for the kinetic temperature 
and the H2 density are ~ 50 to 200 K and 10^-^ to lO'*" cva^, 
respectively. 

The CO abundance per velocity gradi ent relates to th e 
line opacities predicted by the LVG models jWeiB et al.ll2007h . 
Higher values of [CO]/(dv/dr) correspond to higher opacities 
for the high-7 transitions. From Figure |2] we see that all mod- 
els predict optically thick conditions and thermalized emission 
in the /upper < 4 lines, consistent with the line SED following 
and supporting our CO column per velocity gradient selection. 

Although the kinetic temperature is poorly constrained, the 
LVG modelling prediction for the CO 1-0 line intensity does 
not strongly depend on the model chosen: Iqo 10 - 0.27 - 0.4 
Jy km s"'. 
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Fig. 2. CO line SED: integrated line flux, /co, normalized to CO 
2 - 1, vs. rotational quantum number Filled squares mark the 
integrated line fluxes for J100038. Solid lines show LVG models 
A, B and C (see text). Dotted lines show the CO line SEDs for 
NGC253 (Gustenetal. 2006), SMM16359 (IWeifi et al.ll2005al) 
and APM08279 (.WeiB et al.,.2007,). Open triangles show the CO 
SED for the inner Milky Way dFixsen et all 1 19991) . The dashed 
line shows the line flux increasing as v^, which would be ex- 
pected for optically thick LTE conditions. 
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Fig. 4. Spectral energy distribution (SED) for J100038. The pho- 
tometric data points (filled circles) range from radio waves to 
X-rays. The lower purple cross shows the ACS / band flux de- 
rived for the host galaxy while the upper one shows the flux for 
the nucleus. The red and blue cu rves are the averag e radio-loud 
and radio-quie t QSO SED s fromlElvis et al.1 (|1994|) . The dotted 
lines show the lElvis et akl d 19941) models, not corrected for ex- 
tinction, and the solid lines show the SEDs corrected with the 
Calzetti et al. (2000) extinction model. The green line shows the 
SED of the typical starburst galaxy Arp220. The thick gray lines 
show our two-component dust model. 
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Fig. 3. \ogx^ distribution for the single component LVG model 
as a function of T^n and n(H2) with [CO] / (dv/dr) = 10"^ pc (km 
s~') Best fits are provided for solutions with log^^^x^ < 0. 



3.2. CO Size and Mass 

The ratio between the brightness temperatures predicted by the 
LVG models, T^, and the observed line temperatures, T^b, yields 
an estimate of the CO emitting region size. The angular source 
size 9s is related to the beam size 9^ through 



0s = 0b X I 



mb 



-1/2 



(1) 



The source sizes we thereby compute for the different fit 
models (A, B and C) and line brightness temperatures agree well, 
6, « 0.17 - 0.20". At the source redshift (z = L8275), 1" coiTe- 
sponds to 8.44 kpc, so that the estimated source size lies in the 
range of 1 .4 - 1 .7 kpc. 

The CO lu minosity can be deriv ed from the CO line inten- 
sity, following lSolomon et al.l (Il997h as 



Leo = 3.25x10^ 5c 



4s 



K km s 



pc 



(2) 



where 5coAv is the integrated CO line flux in Jy km s Dl 
is the luminosity distance in Mpc and Vobs is the observed line 
frequency in GHz. Our predicted CO flux /i_o = S i-qAV = 
0.27 - 0.4 Jy km s"' yields L\_q = (4.5 - 6.7) x 10^" K km s"' 

2 

pc . 

The CO 1-0 luminosity is commonly used to estimate the 
molecular gas mass. Adopting the conversion factor between CO 
luminosity and H2 mass, a = 0.8 Mf-, (Kkm s~' pc^)~' , that was 
derived for local ULIRGs dDownes & Solomo ni l 19981) . we find 
an H2 mass of (3.6 - 5.4) x lO'" Mq. 

From the size estimate of the dense gas region and the av- 
erage CO line FWHM, Av = 417 km s"', we can estimate a 
dynamical mass contained in the emitting region. For a radius of 
0.75 kpc we infer a dynamical mass Mdyn = 3 X 10"' sin-2(i) 
Mq. For an average inclination angle / = 30°, this leads to 
1.2 X 10"Mq, which agrees well with the CO luminous mass 
estimate. Assuming a disk-like geometry would lower the dy- 
namical mass estimate by a factor 2/7r. 

The molecular gas mass can thus account for a large fraction 
(w 30 - 40%) of the dynamical mass out to 0.75 kpc. This re- 
sult suggests a somewhat larger fraction than in local ULIRGs in 
which the molecular gas ma ss represents 16% of the dy namical 
mass in the nuclear regions ( 'Downes & Solomonlll998l) . 
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Table 2. Infrared and radio fluxes of J100038. 



Band 


Flux density 


Unit 


IRAC 3.6 yum 


116.1 ± 0.2 


fiJy 


IRAC 4.5 fim 


175 ±0.4 


fiJy 


TO A C O . ,™ 

IKAC J. 5 //m 


/JO.l ± 1.1 


fjjy 


IRAC 8.0 fim 


356.4 ± 2.3 


fiiy 


MIPS 24 fim 


1.43 ±0.1 


mjy 


MIPS 70 fim 


7±2 


mjy 


MIPS 160 fim 


50± 15 


mJy 


Bolocam 1.1 mm 


5.6 ± 1.9 


mJy 


MAMBO 1.2 mm 


4.8 ± 1.0 


mJy 


VLA 1.4 GHz 


237 ± 37 


fiZy 



' IRAC and MIPS aperture corrected (total) flux densities are from the 
S-COSMOS data ( S anders et al.. 200 7); Flux densities at 1.1 mm and 
1.2 mm are from Aguirre et al. (2006) and Bertoldi et al. (2007), 
respe ctively; Radio flux is from VLA-COSMOS tSchinnerer et al.i 
I2007h . 



SMGs and high-redshift QSOs have typical H2 masses, 
M(H2) of (3.0 ± 1.2) X 10^° M(?, within a 4 kpc diame- 
ter (iGreve et al l l2005t ISolomon & Vanden Bout 2005), con- 
sistent with the value found in our source. Dynamical 
masses of SMGs and high-redsh i ft QSOs range between 1.0 
and 2.0x10" Mrr, dGreveet alJ 120051: iTacconi et alj l2006t 
ISolomon & Vanden Boutl2005h . which also agrees well with the 
result exposed above. We note that our dynamical mass estimate 
has been computed assuming a radius of 0.75 kpc unlike the case 
of typical SMGs and high-redshift QSOs, for which it is com- 
puted with a radius of 2 kpc. Scaling to this radius, i.e. adopting 
a continuous CO distribution, we find that the dynamical mass 
for J100038 would be 1.6 x 10" Mq, still in agreement with 
values for SMGs and high-redshift QSOs. 



3.3. Dust Continuum 

To study the dust properties of J100038 we use ax^ minimiza- 
tion procedure to fit a 2-component gray-body spectrum to the 
5 IR photometric data points from 24 fim to 1.2 mm (observed 
frame). A single component gray-body spectrum does not pro- 
vide a good approximation to all points. We do not assume op- 
tically thin emission (dust optical depth, Ty <c 1) but use th e 
complete expression for the flux density (see lWeiB et alj|2007b : 



(i+zy 



(3) 



where By(T) is the Planck function, Tjust is the dust temperature 
and Tbg is the cosmic background temperature at the source red- 
shift, Tbc = 2.73 X (1 -Hz) K. We define the apparent solid angle 
subtended by the source as Q = nidolDp}^, with cIq being the 
equivalent source size which we assume to be 1.5 kpc based on 
the derived CO size, and Da the angular distance at the source 
redshift. The dust optical depth is 



Din 



(4) 



where Mdust refers to the dust mass, and the dust absorption co- 
efficient has t he form k(v) = KQjv/vo'f. W e adopt an emissivity 
indexyS = 2.0 (Priddey & McMahon 2001) and kq = 0.4 cm^ g"' 
at 250 GHz (Kruegel & Siebenmorgen 1994). 

The best fitting model leads to a dominant cold component 
with Tc = 42 ± 5 K, and a hot component with Th - 160 + 25 



K. We calculate a dust mass, Mc for the cold component of 
(1 .2 + 0.4) X lO'' Mq; however the mass implied for the hot com- 
ponent is only Mh - (4 + 2) x lO"* Mq. Therefore, the cold 
dust component accounts for almost all the dust mass. We note 
that for a dust mass of ~ 1 ■ 10'' Mq and a size of 1.5 kpc, the 
emission becomes optically thick (tv > 1) at /I < 130 fim (rest 
frame), supporting our model selection. 

The far-lR luminosity can be computed from the modeled 
flux density as 



ipIR — 



Sydv. 



(5) 



Integrating from 50 to 1000 fim dOmont et al.ll200lh . we find 
LpiR - 8.5 x 10'^ Lq. Although the contribution from hot dust is 
only LpiR, H = 5.1 X 10"^ Lq, most of it is provided by the cold 
dust component with Lfir, c = 8.4 x 10'^ Lq. We note that the 
selection of different integration limits imphes different values 
for the luminosity, in particular for the hot dust componen t. If 
we integrate from 8 to 1000 fim dSanders & Mirabe]|ll996 '). we 
obtain Lfir = 1.3 x lO'-' Lq, Lfir,c = 9.4 x lO'^ Lq and Lfir.h = 
3.7 X IO'^Lq for the total, cold and hot component luminosities, 
respectively. 

A substantial difference seems to exis ts in the dust tem- 
peratures observed in SMGs ( Tdust ~ 35 K; iKovacs et al.ll2006l : 
JPope et al. 2006; Blain et al] |2004 and high-redshift QSOs 
(Ldust -45 K: IOmont et al]|2003l:lBeelen et alfcooeHWang et all 
'20081). This difference could possibly be explained by the fact 
that starburst galaxies with warmer dust would have their dust 
emission shifted to shorter wavelengths makin g them more dif - 
ficult to detect at submillimeter wavelengths (Blai n et a n i2004 . 
If this difference is real, it would imply that the dust properties 
of J100038 are more consistent with high-redshift QSOs. 

3.4. Star Formation Rates and Efficiencies 

Assuming that the AGN does not contribute significantly to the 
heating of the cold and hot dust components and that most of 
the dust emission is produced by starburst activity, the far-lR lu- 
minosity implies an estimate for the star formation rate (SFR). 
lOmont et al.l (12001 ) found the following relation between the 
far-IR luminosity and the SFR, 



SFR[M0yr-'] =5 



.^FIR 



10"' Lq 



(6) 



where 5 is a function of the initial mass function (IMF) and vary 
in the range 0.8-3. For 5 = 2, we then find SFR ^ 1700 Mq 
yr \ well in a greement with the value found using the IR based 
estimator froni lKennicut^ d 1 998h . SFR x 2150 Mq yr\ 

The ratio between the far-IR luminosity and the CO lumi- 
nosity can be interpreted as a star formation efliciency, SFE 
= Lfir/l;.q = 130 - 190 Lq (K km s"' pc^)"', or in terms 
of the gas mass implied by the CO luminosity, SFE = 160 - 
235 LqMqI. 

This value is similar to the average value found for loca l 
ULIRGs, (SFE) = (180 + 160) Lq M^ dSolomon et al.lll997h . 
although it appears to be slightly low er than tho se obtained for 
SMGs, (SEE) = 450 + 170 Lq M^ (Greve et al.. 2005) . 

The luminosity ratio can also be interpreted as a gas deple- 
tion time, TsF = M(H2)/SFR ^5x 10'° Mo/1700 Mq yr"' = 30 
Myr. Others (e.g. (irc\c et al. 2005) have estimated the gas de- 
pletion time in SMGs as ~ 16 My r, with a range of 10 - 100 Myr 
dSolomon & Vanden Boutir2005h . well in agreement with what 
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we find for J 10003 8, but shorter than that found in ULIRGs or 
normal spirals. 

3.5. Optical Morphology 
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Fig. 5. Optical morphology of J100038. The HST ACS / band 
image (< 0.1" resolution) is shown with contours scaled in mag- 
nitude. The top panel shows the original image centered at the 
central component C, which coincides with the QSO and X-ray 
source. Two additional components are seen toward the east (E) 
and west (W). The bottom panel shows the nucleus substracted 
HST ACS / image centered at the same position. In addition to 
the E and W components, patchy structure is seen in the central 
component. 



The Hubble Space Telescope Advanced Camera for Surveys 
(HST ACS) / band image of J100038 (Figure [S]) shows a com- 
plex morphology with three discernable components. For the 
bright central component C, Figures |4] and |6] show the SED and 
optical spectrum, respectively. 

The western component is offset by ~ 2" from the central 
source and shows no sign of a connection to it. Although the 



Subaru / band images previously shown bv lBertoldi et alj (l2007h 
are only able to separate C and W, the HST imaging interestingly 
shows a third eastern component (E). 

Since the three components are well aligned, it is suggestive 
that E and W could be gravitationally lensed images. It could 
also be possible th at all three compo nents are lensed images as 
in e.g. APM08279 dWeifi et al.ll2007l) . However, the optical pho- 
tometric redshift of 1 .4 - 1 .5 for W, w hich is resolved in all op- 
tical images (see lBertoldi et al.ll2007l) . is consistent with that of 
C, with no significant secondary solution at higher redshift, thus 
ruling out the possibility that E and W is a background object 
lensed by C. Imaging spectroscopy would much help to estab- 
lish a dynamical relation between the three components, which 
could support a merger history as the cause for the starburst and 
nuclear activity. 

3.6. Spectral Energy Distribution 

Figure|4]shows the J 10003 8 spectral energy distribution from the 
radio regime to X-rays. For comparison, we show the (redshifted 
to z = 1-^3) SEP of the proto-typical starburst galaxy Arp220 
(ISilva et al.ll 19981) . which fits well from radio to IR. The match 
strongly supports the conclusion that most of the radio to far-IR 
emission of J100038 is produced by star formation. 

We also comp are the observ ed photometry with the average 
QSO SEDs of Elvi s et al.l (11994 . While the near to mid-IR pho- 
tometry of J 10003 8 is well matched by both the radio loud and 
radio quiet QSO templates, at higher frequencies we need to ap- 
ply a reddening correction to obtain a good fit. Although the 
radio-loud model would fit the millimeter flux, it is clearly in- 
consistent with the low observed radio flux. To match the model 
at the higher frequenc ies, we applied an e xtinction correction 
based on the model of Calzetti et al. I (12000) for starburst galax- 
ies. A correction with Ay ^ 0.7 yields a good match between the 
models and the UV/optical photometry, and for the radio-loud 
model even for the X-ray flux. Given the redshift (z = 1 .8275) 
and the best fit reddened template, we derive an extinction cor- 
rected absolute magnitude in the Subaru V band. My of -24.2 
magnitudes. 

At near-IR and shorter wavelengths, the photometry cor- 
responds to the central plus eastern component, but at longer 
wavelengths, the western component cannot be distinguished 
any longer In the optical and radio regimes, the flux of com- 
ponent W is only 20% of that of C, but differential extinction 
could in principle be the cause of this difference. 

Based on an empirical point-spread-function created from 
isolated stars in the COSMOS field in the vicinity of J100038, 
Jahnke et al. (in preparation) studied the contribution from the 
AG N and host galaxy t o the optical emission in J100038 (see 
also I Jahnke et al.l 120041) . They find that the host galaxy is re- 
solved and has irregular structure outside the central 0.5", where 
nuclear residuals might dominate, as shown in Figure |5] 

The host galaxy contributes with a ~ (8 + 2)% to the to- 
tal source flux. Therefore, most of the emission at optical wave- 
lengths (92 + 1 % of the total flux) arises from an unresolved point 
source, with /Acs.host = 23.07 + 0.1 magnitudes for the host and 
^Acs.nucieus - 20.43+0.02 magnitudes for the nucleus. This gives 
^Acs.totai - 20.34 + 0.02 magnitudes in total. Component E con- 
tributes a 0.8% to the total nucleus+host flux, /acs.e = 25.6+0.2. 

Figure |4] shows the ACS / band fluxes for the host galaxy 
and nucleus. As mentioned above, the main flux contribution at 
this wavelength is provided by the latter. This is also true for the 
flux in the rest of the optical bands, however it is interesting to 
note that the / band flux density of the host galaxy component 
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Table 3. Summary of derived properties of J100038+020822 



7600 7800 8000 8200 

Observed Wavelength [A] 



Fig. 6. Top: optical Magellan spectrum of J100038 dTrump et alJ 
l2007h . The broad line at A « 7900 A is identified as MglL 
Bottom: close-up of the Mgll line. The smooth curve shows a 
Gaussian fit. The red horizontal error bar represents the position 
of the Mgll line at the average CO redshift (z - 1.8275) and CO 
Une width (~ 400 km s"'). 

matches the optical emission of the Arp220 SED, suggesting that 
it is heavily absorbed by dust and hosts starburst activity. 

3.7. Black Hole Mass Estimate 

To estimate the mass of the central supermassive black hole one 
commonly assumes that the motion of the emitting gas around 
the black hole is virialized (Peterson & Wandel 2000). Then the 
black hole mass is related to the Keplerian velocity of the broad- 
line region (BLR) gas, Vblr, and to the BLR radius, Rblr, 
through Mbh ^blr^blr- The BLR velocity is proportional 
to the line width of the Mgll line (FWHM[MgII]) and the BLR 
radius has been obs erved to correlate with the UV continuum 
emission at 3 000 A (|McLure & Jarvisir2002h . Combining these 
relations (see iMcLure & Dunlopll2004l) . 



M 



ML3ooo \»-^^ / FWHM[MgII] 
Mq ■ ^llO"W/ \ kms-1 



(7) 



where /IL3000 is the luminosi ty at 3000 A (rest frame). In the 
optical spectrum of J100038 dTrump et al.ll2007h . the Mgll Une 
is prominent and allows for an estimate of the black hole mass 
(Figure|6l). Assuming a Gaussian shape for the Mgll doublet we 
find a FWHM of 150.6 + 37.8 A or 5714 + 1442 km s"'. From 
the average radio-quiet QSO SED fit we estimate that the rest 
frame (extinction corrected) 3000 A luminosity is 8.71 x lO-'^ 
W, resulting in a black hole mass of Mbh = (1.7 + O.8)xlO^M0. 

The maximum luminosity that can be reached by an ac- 
creting black hole of this mass, or Eddingto n Umit, is LEdd = 
5.6 X 10'^ Lq. Assuming a typical QSO SED (iElvis et al.lll994 . 
we obtain a bolometric luminosity of Lboi ~ 3 x lO'^ Lq - 
0.54 x LEdd which implies a moderate accretion rate (far from 
the Eddington limit). 

4. Discussion 

4.1. Comparison of Excitation Conditions 

The observations of the CO lines in J100038 show that its molec- 
ular gas is somewhat less excited than observed in local starburst 



Property Name 



Value 



M(H2) 


Molecular Gas Mass (10'" Mq) 
Dynamical Mass (10'" sin"-(0 Mq) 


3.6-5.4 


Mdyn 


3.0 




Dust Mass (10' Mq) 


1.2 


SFR 


Star Formation Rate (10^ Mq yr"') 


1.7 


SFE 


Star Formation Efficiency (Lq Mq) 


190 


TSF 


Gas Depletion Lifetime (Myr) 


30 


My 


Optical Absolute Magnitude 


-24.2 


Mbh 


Black Hole Mass (10' Mq) 


1.7 



or high-redshift QSOs. We find that the turn-over of the CO line 
SED occurs between the CO 5-4 and CO 6-5 transitions 
whereas in most local starbursts/AGNs, or high-redshift QSOs 
studied to date, the peak of the CO line SED is typically lo- 
cated between the CO 6-5 and CO 7-6 transitions. Examples 
of such high er excitation CO e missio n in the local Uni verse 
are NGC253 ( iGusten et al.1 12006'). M8 2 dWeiB et al.l|2005bl). and 



at high-redshift, J1148 +525 1 (Bertold i et alj|2003l: IWalter et all 
2003), BR1202-0725 ("CarilH et al. 2002': ' Riechers et al.ll2006^ 
and APM08279H-5255 ( Weifi et al. 2007 ). The only SMG for 
which the CO SED has been traced over its peak is J 16359+66 12 
dWeiB et al.l l2005ah . which interestingly, shows a similar CO 
excitation to J100038. Other examples of lower excitation are 
found toward th e centers of the l ocal starburst/AGN of Circinus 
and NGC4945 dHitschfeld et al.ir2008,). or in the main starburst 



region of the Antennae galaxy ( Zhu et ai1l2003l) . The lower ex- 
citation of the molecular gas observed in these cases, and in par- 
ticular towards J100038, is likely produced by the relatively low 
H2 density values (n(H2) < 10"* cm-'). Furthermore, the best solu- 
tion derived from the LVG analysis indicates a moderate kinetic 
temperature (~ 95 K), suggesting that the AGN does not con- 
tribute strongly to the heating of the molecular gas. Conversely, 
this may imply that most of the heating is produced by star for- 
mation, as it is also suggested by the similarity of molecular 
gas conditions (Lyn, n(H2)) observed between J16359+6612 and 
J100038. 

4.2. J1 00038+020822: A Starburst-QSO Composite 

The general picture for QSO and stellar spheroid formation is 
based on the merger of two gas-rich disk galaxies. The large 
amounts of gas and dust involved in these mergers provide the 
fuel for infrared luminous starbursts to occur, and the gas inflow 
into the inner regi ons of the galaxy hkely feeds the central mas- 
sive black hole ( Mihos & Hernquis3ll994l: iBarnes & Hernquis3 
1996). Observations of local ULIRGs indicate that the most lu- 
minous phase occurs close to the final stage of t he merger, when 
both galaxy disks overlap dSanders et al.lll988allbHVeilleux et al.l 
1999). 

When the central massive black hole has reached a sufficient 
size and lumin osity, feedback from an active nucleus phase dissi- 
pates the dust dDi Matteo et al.ll2005h and an optically luminous 
QSO emerges. As the gas is being expelled, the feeding of the 
QSO ceases, stopping its activity, and the system relaxes into a 
spheroidal ga laxy hosting a (quie scent) supermassive black hole 
at its center CHopkins et alj|2066l) . 

One important question in galaxy evolution is whether this 
scenario connecting gas-rich starburst galaxies and optically 
bright QSOs constitutes a rule for most of these systems or 
whether it is a phenomenon that applies only to the most lu- 
minous and massive objects. This evolutionary connection has 
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been studied with vast supp orting evidence in the local Universe 
(e.g. ISanders et al.l Il988ah . However, for the crucial epoch of 
star formation and QSO activity (1 < z < 3) only a few stud- 
ies have so far been undertaken. T he study of submillimetre 
jPage et alj2004l:IStevens et al.l2005l) and CO line ("Cop pin et all 
120081) emission of absorbed and unabsorbed high-redshift QSOs 
indicates that the QSO phase could be preceded by a SMG star- 
burst phase and suggests that submillimetre selected QSOs rep- 
resent ideal objects for studying these transitional cases at z > 1. 

J100038 appears to be in a transition between the starburst 
and QSO phases, as a compari son of its propert ies to those 
of a sample of typ ical SMGs (iGreve et alj 12005) and QSOs 
dSolomon & Vanden Bout 2005.) shows. We start by describing 
its QSO properties followed by the features that classify this 
source as a starburst galaxy. 

J100038 shows an optical spectrum typical of a Type-1 
AGN (Figure ^ with a prominent Mgll broad emission line. 
It is relatively luminous at optical wavelengths (My ~ -24 
magnitudes, extinction corrected), shows mild optical extinc- 
tion with Ay ~ 1 and its X-ray emission indicates it is a 
heavily absorbed QSO (logA^(H) = 22 - 23 cm"^). Moreover, 
its optical to mid-IR SED resemble s that of t ypical QSOs 
(Figure |4|, well in agreement with lElvis et al] (11994 ) tem- 
plates. In addition, we derived a diameter of ~ 1.5 kpc for 
the CO line emitting region. This size is consistent with val- 
ues observ ed in high-redshift QSOs which range between 1 
and 3 kpc dWalteret all 120041: ISolomon & Vanden Boutll2005t 
iMaiolino et al .*2007') a nd similar to those found in local UL IRGs 
(~ 0.5 kpc; [bownes & SolomonI [19981 ISoiferet al.ll2000l) . but 
smaller t han the di ameters for the CO emitting region in SMGs 
(< 4 kpc: lTacconi et al. 2006) . 

On the other hand, J100038 shows distinctive features of 
on-going starburst activity. The optical morphology of J100038 
seen in the HST imaging is suggestive of a recent merger event 
(Figure |5]). While most of the emission is concentrated in the 
point-like central source (C), the faint emission from the eastern 
source (E) may be indicative of a tidal tail, hinting at a past in- 
teraction or merger, although it could also imply the influence of 
gravitational lensing. The low number density of HST /-band 
sources in the surroundings of J100038 , p{I < 25.5) - 62 
arcmin"^, implies that the probability of chance association be- 
tween E and C is only 4.1%. If both sources (E and C) are phys- 
ically related, the projected distance between them would be 7.5 
kpc. A large fraction of the optical emission from the host galaxy 
in J100038 could still be absorbed by surrounding dust (Jahnke 
et al, in prep.), and could be hiding the actual link to the eastern 
component. We note that the optical emission of the host galaxy 
is consistent with the obscured SED of the prototypical local 
starburst galaxy Arp220. Indeed, the optical morphology of this 
system is strikingly similar t o that observed in loca l ULIRG/PG- 
QSOs (IRAS 05189-2524: ISanders et all Il988bt ISurace et alJ 
Il998h . Furthermore, the far-IR to radio SED of J 10003 8 agrees 
very well with that of Arp220 (Figure |4|, and does not follow 
the typical behavior of radio-loud or radio-quiet QSOs at these 
wavelengths. In fact, its radio to far-IR spe ctral index, when us ed 
as a redshift indicator results in z = 1.9 ( Bertoldi et alj|20()7l) . a 
strong indication that the far-IR and radio emission are both pro- 
duced by star formation. As mentioned before, the analysis of the 
molecular gas physical conditions suggest that the AGN plays a 
moderate role in the gas heating and hints that the heating may 
actually be dominated by starforming regions. 

All the evidence exposed permits to well fit J 10003 8 in the 
Sanders et al. scenario. Studies of these transitional cases in the 
local Universe have found that about 30 - 50% of the most lu- 



minous ULIRGs (Lp iR > 10'^^ Lq) show broad emission lines 
dVeilleux etal.'1999') and the presence of compact nuclei in 40% 
of all ULIRGs (Scoville et al. 2000; Soifer et al. 2000) as well as 
large molecular gas reservoirs dEvans et al.ll2002ll2005h . even in 
the late stages of this evolutionary sequence. All these properties 
apply to J100038. Furthermore, the presence of both starburst 
and QSO attributes imprinted in the galaxy SEDs have largely 
been observed in nea rby objects, constitut in g the basis of this 
evolutionary scenario dSanders et alJl988allbL[][989l) . These tem- 
plates compare favorably with the SED observed for J100038, 
strongly suggesting that this object is evolving from a starburst 
to a QSO. In fact, it is possible to classify it in a stage between 
the "warm ULIRG" and the "i nfrared excess" QSO phases fol- 
lowing the mentioned scenario dSandersll2004 . 



5. Summary 

We detected CO emission from the millimeter selected AGN 
J100038+020822 at z = 1.8275 in the COSMOS field. 

J100038 is a millimeter-bright, blank-field selected galaxy 
classified via its optical spectra as a Type-1 AGN. The CO line 
intensities peak at the 7 = 5-4 transition and an LVG analysis 
of the CO SED finds that the molecular gas in this source is less 
excited than in typical high-redshift QSOs (e.g. J1148+5251, 
BR1202-0725 and APM08279-H5255), more similai- to what is 
found for SMGs (r^n ~ 95 K, n(H2) = 10^^ cm^^) such as 
J16359-I-6612. The CO emission seems to be concentrated in 
the central kpc, as the comparison between the modeled bright- 
ness and observed main beam temperatures shows. The implied 
molecular gas mass of (3.6 - 5.4) x 10"^ Mq could account for 
a substancial fraction of the implied dynamical mass within this 
radius (< 0.75 kpc). 

Assuming a gray-body far-IR spectrum we derive a dust 
mass of 1.2 ■ 10^ M0 and a dominant cold dust component with 
Tc = 42 + 4 K. The broad Mgll line allowed us to estimate the 
central black hole mass to 1.7 x 10^ Mq. 

Although the molecular gas and dust properties are similar to 
those of typical SMGs and some ULIRGs, the SED from the X- 
rays to the mid-IR is typical of a QSO. The optical morphology 
of J100038-1-020822 is complex, showing evidence of a possible 
interaction or merger Its shared properties of starburst and AGN 
suggest that this composite gala xy fits the evol utionary scenario 
of starburst to QSO proposed bv ISanders et alj dl988a) . 
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